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A GENERAL THEORY OF ASSOCIATION 


BY 
ROBERT GINELL 


Introduction 


When we have a container holding a large number of particles, it has 
long been recognized that the simpler laws of physics and chemistry do 
snot hold exactly if the distance between the particles is too small or, in 
other words, if the particles are too concentrated. The explanation ad- 
vanced is that the residual surface forces present in every particle act 
as an attractive force, compelling some of these particles to agglomerate 
into larger aggregates. While the life span of any one of these aggregates 
may be very short, having such a large number of particles makes the 
number of aggregates present at any instant quite a significant fraction of 
the total number of particles. This fact, of course, would introduce de- 
viations from the theoretical values into our determinations, because 
most of our theories involve, in one way or another, a count of the act- 
ual number of particles present. Qualitatively, this has been recognized 
for many years and, for this reason, determinations are extrapolated to 
infinite dilution or zero pressure, where such association effects, which 
are a function of the concentration, are assumed to disappear. 

- Many attempts have been made to treat this problem quantitatively, 
the commonest being to assume that the deviations are all due to the 
formation of dimers or 2-mers; occasionally to 2-mers and 3-mers or, 
perhaps, even to the presence of 2-mers, 3-mers, and 4-mers.! As each 
additional species was postulated, an additional constant was added. 
Hence the process consisted essentially of seeing how many constants 
were needed to fit a particular variation in some property. In essence, 
this method is always used. If, instead of a purely empirical approach, a 
consistent theory describing the relationships between these camstants 
is offered, then the choice no longer becomes merely an exercise in 
curve fitting. An approach of this type has been to assume that all pos- 
sible sizesof aggregates form, including 2-mers, 3-mers, 4-mers, . 
j-mers; and, in the simplest case, that the rate of formation of any (j+r)- 
mer froma j-mer and an r-mer would proceed with the same rate con- 
stant,? e.g. 


) k ky 
Be + n; 7 ee = K (1) 
F kp ky 
" 5 
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where n, = concentration of particles of degree of association x 
k, = specific rate constant of growth 
k, =specific rate constant of degradation 
K =equilibrium constant 


Here the explicit assumption made is that the same sort of bond is} 
formed or destroyed in every reaction. This assumption probably limitsy 
this case to linear association. An identical approach with, however, the 
the limitation that growth and degradation proceed only through the 
addition or subtraction of l-mers, that is, that the only reaction possible 
is k 


ft PA. 
pet ete; emma UES (2) 


b 


linear association and all species were postulated, but where two cong 
stants were introduced following the reactions 


imine (3j 


IV 
iS) 


pee it eee rt jti j @ 
Calculations using this latter theory have shown that, in certain caseg 
the answers derived give quite satisfying representation of the associatiog 
phenomena.5 The problem involved in using a larger number of constant| 
in the case of linear aggregation has also been discussed but not solvel 
sincea solution involves obtaining as-yet-unknown relationships betwee} 
the multitude of constants which may be postulated. : 

In this paper, the bulk theory of association is presented and a gener 
al solution of the equations at equilibrium is given. The term, bulk theo 
is used to differentiate this theory from the various linear theories, an} 
it covers the linear theory as a special case. 


Theory 


In order to be able to give a suitable picture of the equilibrium 
Situation, we must, however, start with a simpler picture and then go t| 
equilibrium conditions. We shall use the kinetic approach. Let us assum t 
that we start with a very large number of identical particles, which are. 
under such conditions that aggregation proceeds. As a first general a 


addition or breaking off of 1-mers. That is, that the reaction we are cont 
cerned with is 
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ye 


i 
ny + n; a Bs (5) 


Lj+1 
where the specific rate constant of addition of a l-mer to a j-mer is &,, 
while the specific rate constant of the decomposition of a (j+1)-mer to 
‘a j-mer and a l-mer is Lay While this assumption is not altogether 
valid, nevertheless it probably represents the chief process occurring. 
The second assumption which can rather easily be removed, but, which, 
in doing so, restricts the generality of the treatment, is to assume that 
the attractive forces surrounding any particle are spherically symmetric. 
Now when certain conditions, notably energy requirements, are ful- 
filled, two colliding 1-mer particles will combine to form a 2-mer. This 
Q-mer will have acertain life span, and, after some time, it will decompose 
to two l-mers. In each case, this decomposition can happen in only one 
way. The equation would be 


ny * ny n2 (6) 
lo 
and the 2-mer particle can be schematically represented by writing 
ft) SST 
Similarly, when a 2-mer unites with a l-mer, a 3-mer will result. Here, 
however, we have two possibilities; either the 3-mer will be linear (in a 
topological sense) or it will be close packed. This possibility may be 
represented schematically as: 


ny—n;,—Ny or ay — Thy 


ae 


By topological equivalence we mean that we do not differentiate between 
two structures which differ only in the angle of the bonds, but count 
structures only as being nonequivalent, in which a different number of 
bonds are present. These structures, 
p= 1y and Aye pats 
ny 

are equivalent since two bonds are present in each case. This equivalence 
is a simplification of the actual picture, since the force fields extend 
quite far; however, in the case where spherically symmetric force fields 
exist, the forces are extremely short-range in nature and this presentation 
is a fair approximation. If the bond angle becomes too acute, then the 
linear form will be transformed to the cyclic form, and hence we must 
also consider the equilibrium between the two structures. If the assumption 
of symmetric force fields is removed and oriented force fields are sub- 
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stituted, then the situation is simplified, but the treatment is not as — 
general. 
In this case, then, we have the following scheme: 


ng? Get ny n1, 
A 
I I I yu 
k ls k; - 
m =2 m —2 m =1 m = (7) 
k3; jeg) 
hii Sie Se SS SSE n,——nNn, 
1 1 if SS 
L3;, m =3 sys 


In this diagram each multiplicity, m, is also indicated. The multiplicity | 
gives the number of ways in which each step can be taken, and each} 
specific rate constant must be multiplied by its multiplicity when it is | 
used in the rate equations. i 

Similarly, in the case of the formation of the 4-mer, the scheme is | 
as follows: . 


10 Geely 


ny —Ny—)_ ny 


Linear 


Close 


Packed 
A number of interconversions of one form to another are also possible, 
but they have not been indicated in the above diagram. 

As can be seen, this type of scheme becomes very complicated as 
higher j-mers are considered. However, an essential simplification can 
be made which will enable the whole problem to be dealt with mathe- 
matically. This is to say that, whatever the individual structures are or 
no matter how the force fields are arranged, there is for every particular 


over-all step a single constant which, however, may have many parts, 
Peer 
flee 


n+ ny ————s ny+1 (9) 


— 


Lj+y 
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where 
ee Ae Tia, Ik PEF a lTT 
ee te hee + 
; 5K; m; ki m, k te oh (10) 
d 
ae I I ul Teg iil 
a 
bias My hes my, big gee Ty ag Had et it (11) 


Each primed term (Roman numeral) stands for the rate of formation of a 
articular variety of j-mer. 


When these over-all constants are determined, then an investigation 
of their components will elucidate the structures of the various possible 
‘-mers and their equilibria. 


Mathematical Theory 


With this simplification, we can now write the rate equations for the 
formation and disappearance of the various species. 


They are: 
dn 3 > 
Beet SE = iy Kkxyny a Ne tie sonia (12a) 
dt a, x x pe) 
x=1 x=2 


dn nk.on (12h) 
a a Ng ee + tye ton; 


as —nyk3n3 + ny kon2 = lang —13n3 (12c) 


MI 


7 —mgkpny + mp kj-snger + Ly+1 0j+1 ~ Ln; (€2j) 
t 


dny 
dt 


(12m) 


= + nyky-1%m-1 ~ lum 


In EQquaTIoON 12a the rate of change of the number of 1-mers (n,) with 
time (t) is compounded of three terms: the first represents the disappear- 
ance of n, due to combination with each species; the second, the in- 
crease in n, due to the splitting of all the higher terms; the third Huse 
ln., occurs because, when an n, splits, it gives rise to 2 n,’s. 
EQUATION 12b represents the rate of change of the number of 2-mers 
(n,) with time. The factor 1/2 in the second term arises from the fact 
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ee gh ee, Oe 


that two l-mers combine to form each 2-mer. EQUATION 12m is the | 
equation for the largest particle present in the system? This set oft: 
equations obeys the law of conservation of mass which may be stated ash: 
M 
Bice a O 
dt = 
These equations are the rate equations for the process. Our interesty 
is primarily in the situation at equilibrium, which is the case only wherly 
all the left hand terms in EQUATIONS 12 are equated to zero. Useful re 
lationships may be derived from this set of equations in the following) 
manner. 
Let us sum EQuATIONS 12b through 12m. We get 


M 
dn n.k.n 


ee 
dt 7 a long (13) 


x=2 


At equilibrium the individual rates and the sum of the rates is zero 
hence we have 


nykyny 


3 ee aie nl] (14) 


from which 
5 ee (15) 


Similarly summing EQUATIONS 12c through 12m and equating to zero| 
we have | 
k 


ng = L, n yng (16} 


Substituting EQUATION 15 in EQUATION 16, we get 


Similarly we find the general formula 


Ths ei ee ki j (18) 
tS gee 
l, l, l, l 1 
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We now define 


ki _ x 

21, Ky k, : 

k2 [+1 = Ky ps2 (19) 
fe = K» 


defining 
j-l 
nj = nyi TT K ages (20) 
x=0 
From EQUATION 20 we get 
M M j-l 
Paes ipa ashes Ply (21) 
<—1 E fot | x=0 
and 
M M j-l 
> > jnyzi TT Ky (22) 
x=1 ja x=0 


The left-hand terms of EquaTIoNsS 21 and 22 are experimental para- 


meters which can be determined. 
M 


> xn, is the actual count of the number of 1l-mer units in the 


x=1 
system, whether they are present as actual 1-mers or they are combined 


as j-mers. It is the number of formula weights, counting the formula 
weight of the 1-mer as one. The following interrelationship is usually 
necessary to get this number in units conveniently handled. 


M 
> xn, 
M x= = Ww (23) 
x = a 
ee hd ala 


<—1 


where now SxW_ is the count of the number of 1-mer units in multiples of 
Avogadro’s avert No and is equal to the weight of sample w, in 
grams, divided by W,, the weight of an Avogadro number of 1-mers. This 
number, which is the stoichiometric concentration, can be called the 
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M 
formolity (moles according to formula). The }>p is the actual count B 
messy 


of the number of particles present. Similarly 


gee (24) 


x=] N 


Saviecn be called the avmolity (Avogadro moles). It can be obtained § 
x=1 \ 
from a measurement on the system; for instance, from the colligative ¥ 
properties. 

From EQUATIONS 21 and 22, we must now determine the value of f 


n, and the values of the various It’s. Let us proceed stepwise. 
Case where only 1-mers are present, This is obvious in its solution, 


although trivial, as few if any systems contain only l-mers except in 
dilute systems. Here 


1 if i 
SS ee De xn, =n, (25) | 


x=] x=1 


and, of course, there are no higher terms. 


Case where only 1-mers and 2-mers are present. Here we have: | 


xn =n +2K n 2 
= 1 2k in, 
Again the solution is obvious and simple since 


- 2 
2 a me 
ny xn, =n, 
1 1 


and 
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This case may be met in certain peculiar situations (usually as a first 
proximation), where we consider particles having only one point of 
ossible attachment. 5 

Case where only 1-mers, 2-mers, and 3-mers are present. The equa- 


tions are: 
3 
> = ; 2 py... 3 
been thon + 
= x 1 fut K,h,n, 
D ce® | (28) 
3 
= 7 2 , 3 
xn =n_ + 2h + 
> r. 7 2 i 3K k,n, 
x=I 


Here we have three unknowns, n,, Kz, K2; hence we cannot determine 
the values from these two equations alone. We must take the values at 
two different concentrations; getting thereby: 


(Zn) =(2,),+ K,@),’ +K,K(n,),? 


4 7 oa , 
(= xn), =(n,), + 2K (n,), + 3K K(n,),? 
(29) 


2 n), = (n,), ee K (a,),’ de KK (a),° 
(X xn,), =(n,), + 2K (n,)," * 3K K,K,(0,).° 


This calculation gives four unknowns, (a) (tt es I ae and, having 
four equations, the set may be solved analytically. 

Cases where more than three species are present. For such higher 
cases, an algebraic solution is always possible. For four species we 
should use three points and have six equations and six unknowns; for 
five species, four points, and have eight equations and eight unknowns, 
etc. However, the algebraic manipulation involved becomes forbidding. 
A simpler approach to this problem is to use matrix methods, which by- 
pass the complicated algebra and resolve the problem into one of simple 


straightforward numerical computation. 


Matrix Solution 


Let us assume that our basic series, having j terms, are 


a) 2: “4 2 7-7” 3 ro “4 Jj 
Ln, =n, +K,a, +h kn, #0, +h ky. Kea, (30a) 
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and 


7 %) “T7 7 j 0 . 
Kon ie gk hoe Eee (30b)e 


7 2) 7 
Xxn =n, +2K,n, +3K Kn, 


For convenience we shall change our notation to 


U = X + aox? + a3x3 + agx4 +... 4ajxi (31)E 
V=xX+ Cox? +.03x3 + cyx4 + ...40;xI (3206 
where 
u; = (Xn,); xX=Nny a, =K, Cy = 2Ky 


XxX, =(N1); a; =K,K, e5e JK 


v; = (Zxn,); Xj=(n,); a; = K,K2_K,_, c, = jK,K, ouelyly 


We have now determined for the j points the values of u; and v,;. Let us|: 

assume that we know the values 
X77 X2,0X35) X4, CLC. 

(The determination of Xi» Xp etc. will be considered in the next section. Jif 

Now for j points, the following set of equations can be written: 

ul 


Xy + AX 17+ a3xz3 +... + ajXyI 


u, Sy gk et Aka? ee Eas (33) 


. 


(= 2 2 4S} j 
u; Xj + AX; + 3X; Fees $aAjzXi! 


In matrix notation this calculation becomes: 


uy = Xp Xow Xk? eae Z 
U2 X92 Soe Koo vee Xoj a2 
u3 X3 xe x33 S08 x3/ a3 
: (34) 
u; 2 3 j 
J x; X 4 Scent A; 
nit j ia ay 


— a 


Since [Xx,i] is a matrix of numbers which is nondegenerate, the equation : 
is formally solved by writing 


Xpe XK) Wek y eee! —1 uy a 1 
X9 X92 X55 Voom! u> a, 
Xo Xg2 Kao Tuy, ees u3 a; 
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The explicit solution which is desired can be obtained by numerically 
computing the reciprocal of the [X;i|] matrix. The product of the first row 
of the reciprocal [X;,i] matrix multiplied by the [u] column matrix will 
be equal to one. The product of the second row of the reciprocal Xi 
matrix multiplied by the [uw] column matrix will be equal to a,, which is 
K,. The product of the third row multiplied by the [u] matrix will be 
equal to a3 which is K ,2 ; from the determined value of K, the value of 
K2 can now be obtained directly. Similarly the rest of the K’s can be 
gotten. There are several well-known methods of reciprocating matrices. ° 
These methods will not be discussed here, but the problem is one of 
straightforward computation. 


The Determination of ny 


There remains the problem of determining the value of the x,’s be- 
fore this solution can be applied. Our original series is 
x teen + a,x! (31) 
This is a finite power series. Now there exists an infinite series such 
that 


u=x+a,x2+a,x3+a 


x =u + bou? + 63034 bgut +... (36) 
It has been shown that such an inverted or reverted series exists even in 
the case when the original series is an infinite power series.’ The co- 
efficients in the reverted series can be determined by substituting EQuUA- 
TIon 36 in EQUATION 31, expanding, collecting terms, and equating the 
coefficients of like powers of u . This can be illustrated for the series: 
u = x + 42x ?4+ a3x? (37) 
x=u+bju2 +b,u3 + byut+bsurt... (38) 


Substituting EQUATION 38 in EQUATION 37 we have 


u=(u+by 1?+b 947+... )+a,(utb2u7+b u%+...)?4a,(atbau +b 37 +...)” (39) 
Expanding, collecting terms, and equating like powers-of u, we have 
u=u 
u?(b,+a,) = 0 
u3( b3+2a2b2+a3) =0 (40) 


u4* [b,ta,(2b,+b, 7) + a; (3b,)] =0 
u5 [bstay(2b,4+2b ab 3) + a3 (32, 495.7 1 0 


etc. 
Since the u’s are not equal to zero, we obtain 
bg =— &2 
b3 = 2anb 2—-a3 = 2a7?—a3 (41) 
bg=- bag? ai 5azaz3 


etc. 
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The coefficients of the reversion or inversion of an infinite power serie sti 
have been determined to the 13th coefficient.* The first few coefficients 


of the reverted infinite series are: 


b2=—€a2 
b3 =+ 2a 27-a3 } 
bg = —5a,3+5a7a3—a4 (42) 
etc. 


As can be seen by comparing EQUATION 41 with EQuaTion 42 thel 
solutions for the first two coefficients are identical. The third terms, 
ba, differ by the factor, ag. This factor, ag, does not appear in EQUA- 
TION 37 (i.e., it is equal to zero), hence it does not appear in the com: 
ponents of the 5 coefficients. If additional terms are added to EQqua-—jJ 
tion 37, then more coefficients in the inverted finite series would be | 
identical to the coefficients in the inverted infinite series. By applying? 
this principle we can derive the coefficients up to the 13th term 
by cancelling the appropriate components in the known coefficients of 
the reverted infinite series. If necessary, further coefficients can be com-p 
puted, although the computation is tedious. Fortunately not more than 
the first few will probably be required. We now have the following set of | 
inverted infinite series: i 


Xj = u,+b 2u; 245 3u 7 34+ Ge, 


X2q = Up+b 2 27+b3u 25+... 


xX; =Uj+bou,;?2+b3uj;5+... 


which are obtained from EQuaTions 33, and the set of inverted series) 
obtained from EQUATION 32. 
Xy => Vytdov, 2+d3v,3+ see 


X2 = Votdov2? 4+d3v254 ... 
(44)! 
=) °F Vj+d2Vv; 2+d3v ;3+ Sac 
where 
d, =—Co9 
pry P34 
d, = P2cs == Cig 
d4 = —5c» 345 C9C3—-Cy4 
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he matrix equation for these two sets are: 


x1 | = [ a7? ay?) 4 a 


X2 et os a ite. bo 
xs t, w,7 a, ... | b, | 
bg 
. ee (45) 


d 

xX = Vy Wy? V,3 1 
X2 V2 V2? v29 d2 
X3 V3 v3? v3? d3 

; d4 

(46) 
: 2 vy.3 

ga Vv; Vv; Vv; — d; 


these [u] and [vy] matrices cannot be reciprocated since they are 
ectangular matrices, with infinite rows. However, since these series 
power series and are convergent, we can use an approximation method. 


The convergence of the u and v series is an interesting problem. 
If the formolity for instance is above one, then the series 


Vy v7", vz? wee 

is of course divergent. To convert this series to a convergent series, 
we must express the formolity and the avmolity in such units, so that 
the largest determined point is below one. Probably it would even be 
advisable to have the values lie well below one, so that the series 
converge rapidly. After the x’s have been evaluated, the formolity 
can again be expressed in the original units. 
Let us assume that x can be represented to the desired precision by 
three terms of the u, or v series; then, taking the values for three de- 
erminations, we have the set of equations: 

X1,= u +b ou 17+b3u 17 

Xg = Uztb ott 2?+b zu 27 (47) 

x3= u3+b 2u37+b3u37 
and 
x1 = vitdov17?+d3v17 
x2 = votd2v27+d3v27 (48) 
x3 = v3tdoVv37+d3v3" 
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In matrix form these sets become: 


X2 | u2 
| x3 | ee 
and _ = 
| at | Fi Wi 

269) V2 

X3 | WG 


uy," 


U2 


u3? 


yo | 1 ] 
use | 22 
u3? | bs! 
v1? ie 
Vas we | 
vat | [aa 


These equations can be solved directly to give: 


— 


fe hee ieee x1 | = ‘iy 
U2 U2? U2 | X2 bo (Sa 

| u3 u3? u3 | cies ft | 53 | 

and Zz a 

| va Vi" vio |4 leexa = | 2 q 

| V2 V2? V2? X2 do (52) 
a V3 us” X3 | ds | 


The reciprocals are easily obtained, since the lu] and [vy] matrices 
consist of numbers, and we get symbolically: 


From EQUATION 51 


r,s; (tj x, Ce cen) Oy 
roe sy 3 X2 bo (53k 
yee NER Lae x3 | nM b3 | 
and from EQUATION 52 
[tiem Pe pin ye fe ang ee 
ly iY9 po X2 d2 (54)' 
Ls l3 m3 P3 X3 | d3 


where r, s, t, and 1, m, p are numerical coefficients. 
We now have from EQUATION 53 


Ty Xy+S 1X ptt 1x3 = 1 (35) 

tf X1+S2xX2+tax3 = bo 
and from EQUATION 54 

lyx1+m1X2+P1x3= 1 


(57) 
Lox y+ 2X t+ Pox3 = dz (58) 
Since b, = — Qo 
and d, = — Cy, = — 2a, 
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e get from EQUATIONS 56 and 58 
(lz-2r) XJ + (m2—2s 2) Xt (p2—2t 2) X3 =0 (59) 


ZouaTions 55, 57,and 59 in matrix form are: 


1 Si ( i xy = ; 7 
ly m1 Py > 7 (60) 
(lz—2r2)(m2—2s2)(p2—2t2)|_ | Xs | 0 
hose solution is 
=| = t1 S1 ty toltl 
X2 ly my P1 1 (61) 
X3 (1—2r2)(m2—2s2) (p2—2t 2) | 9 | 


The explicit solution for x,, x2, and x3 follows directly, once the 
xeciprocal has been determined. 


We now must repeat this computation, using 4 or 5 terms and, cor- 
respondingly, 4 or 5 experimental points, so as to insure that we are 
getting the x’s to the desired degree of precision. Once we have de- 
ttermined how large a matrix is necessary in order to get the desired 
significant figures in x, then, using the same order matrix, we can repeat 
the computation of the rest of the j points so that we obtain all the j 
values of x. 


Now, since we know the values of x, we can solve either these 
matrices forthe first few K coefficients, or, what is probably much simpler, 
use EQUATION 35 to obtain the values of the various K’s. 


Discussion 


The great advantage of this treatment lies in the fact that it isa 
method whereby the equilibrium constants of association may be de- 
termined empirically from the experimental parameters with few if any 
assumptions, explicit or implicit. Explicitly, we made one assumption 
in setting up the rate equations, that is, that growth and degradation 
proceed through the mechanism of 1-mer addition or breaking off. However, 
this assumption has no bearing on the situation as it exists at equilib- 
rium. It has been shown? that no matter what mechanism of growth and 


decay is postulated, at infinite time or equilibrium, the degree of 
e constants. Since the degree of 


Sn, is also depend- 
on independent 


association is a function only of the rat 
association is defined as equal to «"x/Sn,, then 
ent only on the rate constants, since > «7, is by definiti 
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of the mechanism. The other explicit assumption, namely, that the force mi 
surrounding the particles were spherically symmetric, was only pertinent: 
where we considered the various forms that could exist of each j-me i 
One implicit assumption or restriction which does nevertheless limit the : 
treatment, is that we assume that no dissociation of 1l-mers occurs. Thi! 
assumption means that we must choose the l-mer very carefully, ang 
that it must be the smallest aggregate which is associating. If this 1-me! 
also dissociates partly, into dissimilar smaller units then the treatmen’ 
becomes in general much more complicated. Nonetheless certain particule 
cases of this sort may be treated. !° I 


i 


J 


the question of what the component rates and multiplicities of the indil 
vidual types are becomes more amenable to investigation. | 


Work is continuing on this problem. 


Summary 


(1) The methods of approach to a bulk theory of association are des) 
cribed and discussed. 


(2) The rate equations for growth and degradation of j-mers, postulat} 
ing that all species can be present and assuming a mechanism of 1-mea 
addition and degradation, were derived. t 


(3) The solution of this set of equations at equilibrium is given) 
using matrix methods. This solution permits the values of the equilibriu 
association constants to be derived empiricaily, free of any assumptions, 
with, however, the implicit restriction that dissociation of 1-mers does! 


not occur. ! 
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